TITLE: 
TYPE: 
CAPACITY: 
DURATION: 



ORDER OF TAPES: 



DATA TAPE: 



PURPOSE: 



SUM CHECK: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER LABORATORY 
STATISTICAL LIBRARY 

KSL 2.50 - 2^8 

Fisher's Z Transformation: Z=| ln(l+r)/(l-r) (SADOI Only) 

Entire Program 

No restrictions 

About to) ms for each three-place number to about k^O ms 

for each twelve -place number 

1. Master tape: stops on 2kO06 

2. Data tape: with black switch at obey, data tape will 
stop at 2kO06 when an "N" is read; with black switch at 
disable, data tape will stop at OF when a J is read. 

3. Additional data tapes may be read by raising the 
white switch. 

The data tape consists of a set of signed fractions scaled 
by 10" 1 and terminated by an N or an NJ. Tfcere must be at 
least one N at the end of the tape. If there are several 
N's separating groups, then it is desirable *> punch a final 
J character so the problem may be run on stop disable. 
The typical data tape is the output of K-8, K-9, or any 
other correlation routine (either a triangul^ or square 
matrix) terminated by an N or NJ. 

The sampling distribution of the correlation coefficient, 
r depends upon the size of the population value, r , wM 
usually is unknown, and also it depends upon the sample 
size, N, Fisher's Z, 

Z = tanh" 1 r = I In (1 + r)/(l - r) 

is so nearly normal for any value of r that it may be 

treated as such for all practical purposes . 

The standard error of Z may be estimated by the formula, 



', which 



S = 
z 



/" 



N - 3 

The chieT uses of Z are to be found in problems of 
averaging coefficients of correlation and in testing 
the significance of differences between r's. 
The program tape is sum-checked during input* A failure 
to input the program correctly will result in an FF stop 
from location 1L5. 

For reference, see: R.A. Fisher, Statistical Methods for 
Research Workers , 10th Edition, Oliver and Boyd, Edinburgh, 
19^8. 



Fisher's Z 



- 2 - 



NOTE: 



The value of Z can range from positive to negative infinity, 
A test is made in the routine such that whenever the 
absolute value of r attains or exceeds .999999 x 10"^, 
then a dash is punched on the tape, and the routine goes 
on to the next r value. 



DAT E December 8, 1958 
SUBMITTED BY w X fft^UP^^ 
APPROVED BY ^ 
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ORDER 



' MOTES 
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FtSHER'S TRANSFORMATION 
Z =i" MT (1+R)/(1-R) 

R Z 

006K 

L5(-) 0020F 
463L 50F 
L5(E1) 4o(e) 
50F 503L 
26(N12) 40F 
L3F 326L 
2213L L521(N1?) 
L0(13) 1020F 
40(E) LO(-) 
40F L3F 
3211L 5010L 

26 (Z) 92131F 

92770F 92131F 
24L 92131F f 
92834F 92135F 
OFF 26L 

OOK 

OOF 00210F 

OOF OOF 

OOF 001023F 

OOF OOF 

OOF OOF 

OOF OOF 

40102 3F L52l(N12) 

OOF 0&1900 0000 OOGQJ 

OOF 000999 9990 OOOOJ 



Title 



14a jj! routine 



Set store address 



KSL 2.50 



Read in data; stop on N or store at 1023 
Test for final J 



Set final atore address at (E) 



24006: stop at N 



OF: stop at final Jj vhite switch for 
new problem 



Initial store of r's 
Final store of r 

Tally counter 
Temporary storage for r 
Temporary storage for In (l-r) x 2 
Test constant for (N12) 

10- 1 

Largest r that will be transformed 



LOCATION 


T" 


ORDER 
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KSL 2.50 


fcbs. 


Rel. 


Sym. 


OOF OOF 






31 




(10) 




Temporary Storage 




32 




(11) 


OOF OOF 




Temporary Storage 




33 




(12) 


OOF OOF 




Temporary Storage 




3^ 




(13) 


kOF L52l(N12) 




Test constant 




35 




(HO 


26 (N) OOF 

OOK 




Transfer order for (N12) 




36 





(z) 


K5F ^221L 




Subroutine to calculate and print Z's 




1 




L5(-) ^23L 










2 




L5(E) *£22L 




Set initial address for r 






3 




1*1 (c) L5F 




Set final test constant 






1* 




1*0 (R) 50F 










5 




5^103F 505L 










6 




26(Pl6) 9296IF 


Print r 






7 




L5(2) 12 (R) 




Test r for maximum size 






8 




329L 92T1^F 


1 If too large, print a dash 






9 




2218L L5(l) 








10 




lo(r) 5010L 




pr 






11 




26(S5)%)(LN) 




In (1-r) x 2~? at (LN) 






12 




L5(l) iMr). 










13 




50F 5013L 




In (l+r)/(l-r) x 2 -5 at N(o) 






1* 




26 (S5) LO(LN) 








15 




1*0F 50(1) 




|ln (l+r)/(l-r) x 10" 1 






16 




75F OO^F 








IT 




5^105F 5017L 










18 




26(Pl6) 92131F 


Print Z 






19 




92515F F53L 










20 




1*2 3L L022L 










21 




323L 22F 










22 




N1(C) L5F 








59 




(S5) 


OOK 




Natural logarithm 




95 




(Pl6) OOK 




Infraprint 




151 




(N12) OOK 
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Abs. 
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1 
2 



Sym, 



(N) 



ORDER 



OOK 

IMN12) 4621 (M12) 
[from (N12) 

Ll(T) L42l(N12) 

563L 2623(N12) 

L5F 4o(lO) 

L51F 4o£.l) 

L52F 4o(l2) 

50F 506L 

26(Z) L 5 |(-) 

0020F 4621 (N12) 

L5(l0) 40F 

L5(ll) 401F 

L5(l2) 402F 

r 
264 (N12) OOF 

00500K 

15(14) 4022(KL2) 

26999f OOF 

00500K 

L3F 346F 
FFF 266F 
KL3579F 282430F 
26L 26lN 
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■*■ 



Appendage for (N12) 



Test for store order in (N12) greater 
than 1023 



KSL 2.50 



Transfer to subroutine to calculate Z's 



Return to (N12) to continue inputting r's 
Interlude to set transfer in (N12) 



Sum check 



